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\Q . Abstract 



The cosmological stochastic gravitational-wave background produced by the mildly 
non-linear evolution of density fluctuations is analyzed, in the frame of an Einstein-de 
Sitter model, by means of a fully relativistic perturbation expansion up to second or- 
Qs, . der. The form of these gravitational-instability-induced gravitational waves is gauge- 

dependent. In the synchronous gauge, where the second-order expansion is most easily 
carried out, the transverse and traceless tensor modes which are produced also contain 
a Newtonian and post-Newtonian piece, whose interpretation as gravitational waves is 
non-trivial. A more direct physical understanding of this background is obtained in the 
so-called Poisson gauge, where it is seen to consist of a constant term plus an oscillating 
c/3 . piece whose amplitude decays inside the Hubble radius. 

a 

> 

X ' 1 Introduction 



In these notes we will review some recent results concerning the generation of gravitational 
waves from the gravitational instability of density perturbations in a cosmological framework. 
The mixing of different modes (scalars, vectors and tensors) is a consequence of non-linearity, 
which already arises in second-order perturbation theory. To this aim we will first consider the 
exact relativistic dynamics of irrotational dust in the Einstein-de Sitter background (i.e. spa- 
tially flat Robertson- Walker geometry, with vanishing cosmological constant), in a synchronous 
and comoving frame. We will then expand the exact equations in this gauge, up to second order 
in the amplitude of deviations from the background solution. A second-order gauge transfor- 
mation is then performed in order to obtain the metric perturbations in the so-called Poisson 
gauge, where the tensor modes are most easily interpreted as gravitational radiation. 

Throughout these notes we use signature +2 for the metric; Greek indices take values from 
to 3, and Latin ones from 1 to 3. 

The line-element will be written in the form 

ds 2 = a 2 (T)g IIl/ dx' 1 dx u , (1) 



where a(r) is the scale-factor of the background flat Robertson- Walker universe and r is the 
conformal time. We deal with a pressureless fluid, with vanishing cosmological constant, so 



the background Einstein-de Sitter model has a(r) oc r 2 . The components of the perturbed 
conformal metric tensor can be written as 
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where xt^ = 0- Indices are raised and lowered by the Kronecker symbols 5 U and 5^, respec- 
tively. The functions ip^ r \ c<4 , (f)( r \ and xjj represent the r-th order perturbations of the 
metric. 



uric. 

We can split the perturbations into scalar (i.e. longitudinal), vector (i.e. transverse) and 
tensor (i.e. transverse and traceless) modes. In particular, uf can be decomposed as 



M <-> frill i (r)± 



(5) 



where cu- r ' ) " L is the vector component, i.e., d l u>^ ± = 0. Similarly, the traceless part of the spatial 
metric can be decomposed as 

= D^ (r)l1 + d a { ; ] x + d jX ? )L + xl; )T , (6) 

where x^" is a scalar potential, xt^ are vector modes, and xtj T tensor modes {d l Xij^ T = 0); 
in the equation above we introduced the symbol 



B il =d i d 1 ~-5 l1 V 2 . 



It is also useful to introduce similar expansions for the density contrast 
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with g(x., r) the matter density and Qbij) = 3/27rGa 2 (r)r 2 its background mean value, and the 
fluid four-velocity tt M , 

1 / 00 1 \ 
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The four-vector it M is subject to the normalization condition u^u^ = — 1; therefore at any order 
the time component is related to the perturbation ip(r). For the first and second-order 
perturbations one obtains [|TJ 



v °(i) = -V'(i) 

,,0 



The velocity perturbation i>K can also be split into scalar and vector parts 

v\ r) = &v\ r) + v\ r)± . 
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2 Relativistic dynamics of irrotational dust in Lagrangian coordi- 
nates 



We start by writing the Einstein's equations for a perfect fluid of irrotational dust in the 
synchronous (g o = — 1> 9oi — 0) an d comoving (u M = Sq/o) gauge. The formalism outlined in 
this section is discussed in greater detail in Ref. 0. 
The line-element is now written in the form 

ds 2 = a 2 (r)[ — dr 2 + ^(x, r)dx l dx s \ , (13) 

with the spatial coordinates x representing Lagrangian coordinates for the fluid elements. 

By subtracting the isotropic Hubble-flow, one introduces the extrinsic curvature of constant 
t hypersurfaces, 

with a prime denoting differentiation with respect to the conformal time r. 

One can then write the Einstein's equations in a cosmologically convenient form. The energy 
constraint reads 

O OA 

$ 2 - tf^i + ^ + n = —6, (is) 

where lZ l j(g) is the intrinsic curvature of constant time hypersurfaces, i.e. the conformal Ricci 
curvature of the three-space with metric and 1Z = W ^ 
The momentum constraint reads 

= *4 . (16) 

where the vertical bar indicates a covariant derivative in the three-space with metric g^. 

Finally, after replacing the density from the energy constraint and subtracting the back- 
ground contribution, the evolution equation for the extrinsic curvature reads 

+ + + \ (f<d* k - + TV 3 - \lZ5^ = . (17) 

Also useful is the Raychaudhuri equation for the evolution of the peculiar volume expansion 
scalar namely 

+ ^^ + —5=0. (18) 

An advantage of this gauge is that there are only geometric quantities in the equations, namely 
the spatial metric tensor with its time and space derivatives. The only remaining variable, the 
density contrast, can indeed be rewritten in terms of gr^, by solving the continuity equation. 
We obtain 

«5(x,r) = (1 + 5o(x))[s(x,t)/s (x)]- 1/2 - 1 , (19) 

with g = det gy. 



3 Perturbative approach in the synchronous gauge 

In this section we will obtain the perturbations of the Einstein-de Sitter cosmological model up 
to second order in the synchronous and comoving gauge. Different approaches to this problem 
have been used so far. The first solution of the second-order relativistic equations has been 
obtained, in this gauge, in a pioneering work by Tomita ||. Matarrese, Pantano & Saez 0, || 



obtained the leading order terms of the expansion, using a different method, based on the so- 
called fluid-flow approach (e.g. Ref. ||). Salopek, Stewart and Croudace JF| used a gradient 
expansion technique to obtain second-order metric perturbations; an intrinsic limitation of their 
method is, however, that non-local terms, such as the non-linear tensor modes, are lost. Russ 
et al. || recently rederived the metric perturbations to second order in the synchronous gauge, 
using a tetrad formalism. 

We start by expanding the covariant conformal spatial metric tensor up to second order in 
the form 

9ij = + $ + \gt? . (20) 
The contravariant metric takes the form 

g * = $s_ g <m_ i g m + g m*gM. (2 i) 



The extrinsic curvature tensor ■ up to second order reads 

= + ^ < f-9 (1)i VS ) '). (22) 

The square root of the metric determinant is 

9 112 = i + \ 9 { f + \g { ? + l(g { ?) 2 - \ 9 ™$ , (23) 

with inverse 

9- 1/2 = 1 " \ 9 { f ~ -/? + \{9 { ff + \ 9 ™$ , (24) 



so that the density contrast reads 
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having assumed as initial conditions t/oij — an d ^o^ = (i- e - = ^cj 1 '*)- 
The Christoffel symbols up to second order read 

T --(a {1)i +a {1)i -a (1) A+-(a (2)i + a {2)l - a (2) A - -a (1)l ( a (1)e + a {1)£ - a {1) A (26) 
1 jk — 2 1 9 j ,k + y fe j yjit J + ^ 1 9 j ,k + y & j yjfc J 2^ e [9 j ,k^ 9 k j 9jk J > \ zo ) 

from which, after a lenghty but straightforward calculation, the conformal Ricci tensor of the 
spatial hypersurfaces is obtained 
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its trace reads 



(27) 
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3.1 First-order perturbations 

We are now ready to deal with the equations above at the linear level. Let us then write the 
conformal spatial metric tensor in the form 

g tJ = 5 tJ +g%, (29) 

where, from now on, we use the subscript S to label synchronous-gauge perturbations. 
We then write 

,g> = -24% + d,^'" 1 + daif + a«r + xg )T , (30) 

with 

,(1)-L (1)T» _ ™ fl)T 



^sr = ^r = ^> = 0, (si) 

indices being raised by <5 lJ . Remind that at first order the tensor modes xf^ are gauge- 
invariant. 

As well known, in linear theory, scalar, vector and tensor modes are independent. The 
equation of motion for the tensor modes is obtained by linearizing the traceless part of the i? 1 ■ 
evolution equation. One has 

xg^ + ixg^-Vxg^-O, (32) 

which is the equation for the free propagation of gravitational waves in the Einstein-de Sitter 
universe. The general solution of this equation is 

x5 1} (x,r) = JLy rf 3 kexp( ? k-x) X «(k,r) e ^.(k), (33) 

where e^-(k) is the polarization tensor, with a ranging over the polarization components +, x, 
and (k, r) are the amplitudes of the two polarizationm states, whose time evolution can be 
represented as 

X^\k,r) = A(k)a^)(^^y (34) 

with ji the spherical Bessel function of order 1 and a CT (k) a zero mean random variable with 
auto-correlation function 

Mk)cv(k')) = (27r) 3 AT 3 5 3 (k + k')<W • (35) 

The spectrum of the gravitational-wave background depends on the processes by which it was 
generated, and for example in most inflationary models, A(k) is nearly scale invariant and 
proportional to the Hubble constant during inflation. 

In the irrotational case the linear vector perturbations represent gauge modes which can be 
set to zero: xf^ — 0. 

The two scalar modes are linked together via the momentum constraint, leading to the 
condition 

4> + 1 V\<"« = 4 1 ' + Iv 2 xS" ■ (36) 
The energy constraint gives 



^ , ' + 40i" , -xS , )+2«g + iv^> 



, (37) 



r 2 



having consistently assumed 6q <C 1. 

The trace of the evolution equation also gives an equation for the scalar modes, 

x g»* + i x g)l' + iv» X g» = -2*if . (38) 

An equation only for the scalar mode Xs can De obtained by combining together the evolution 
equation and the energy constraint, 



V 2 



\s . (39) 



r 2 



On the other hand, by linearizing the solution of the continuity equation, we obtain 

J?> = *,-iv»(xg» -,&"), (40) 

which replaced in the previous expression gives 

4 1) " + ^4 1)/ -4 7 rGa 2 ^4 1) = 0. (41) 

This is the equation for linear density fluctuations (e.g. Ref. [[J) whose general solution is 
straightforward to obtain. 

The equations above have been obtained in whole generality; one could have used instead 
the well-known residual gauge ambiguity of the synchronous coordinates (e.g. Refs. 0, ||) to 
simplify their form. For instance, one could fix Xq 1 '"' so that V 2 Xsd" — — 2<5 , so that the Xs 
evolution equation takes the same form as that for 8g. With such a gauge fixing one obtains 

X^'W) =X+(x)r 2 + X _(x)r- 3 , (42) 

where x± se ^ the amplitudes of the growing (+) and decaying (— ) modes. In what follows, we 
shall restrict ourselves to the growing mode. The effect of the decaying mode on second-order 
perturbations has been considered by Tomita || and by Russ et al. |§ and will not be studied 
here. The amplitude of the growing mode is related to the initial peculiar gravitational potential, 
through x+ = — Iv 9 ; where in turn, is related to 5 through the cosmological Poisson equation 

W(x) = JU(x) . (43) 



Therefore, 



D^'^-y^-^VV) . (44) 



The remaining scalar mode reads 

4 1) (x,r) = ^(x) + ^VV(x). (45) 

In what follows we will also drop the tensor mode Xi^ T as a possible source for second-order 
perturbations, i.e. we will restrict ourselves to initial growing-mode scalar perturbations. The 
effect of these modes on second-order perturbations has been studied by Tomita JTOl O and 
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3.2 Second-order perturbations 
Let us start by the definition 

9 { sl = -2</4% + ^ , (46) 

with xs )l j ; = 0. 

The technique of second-order perturbation theory is straightforward: with the help of 
the relations reported at the beginning of Section 3 [Eqs.(20) - (28)], we first substitute the 
expansion above in our exact fluid-dynamical equations (momentum and energy constraints 
plus evolution and Raychaudhuri equations) obtaining equations for g^, with source terms 
containing quadratic combinations of plus a few more terms involving Sq. Next, we have to 

solve these equations for the modes (f)^ and x^slj m terms of the initial peculiar gravitational 
potential ip. 

Let us now give the equations which govern the evolution of the second-order metric per- 
turbations. 

The second-order Raychaudhuri equation reads 
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The second-order energy constraint reads 
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The second-order momentum constraint reads 
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(48) 



z< Ps,j + -jXsj ,i - 9s [9sjk,i - 9sik,j J + 9s ,i9sjk - Tfls ,j9sik ~ ^si ,k9sj ■ l 4y J 
Finally, the second-order evolution equation reads 
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+ ^S fe ' m ^V (50) 

The next step is to solve these equations. In these calculations, we can make the simplifying 
assumption that the initial conditions are taken at conformal time r = 0. One can start from 
the Raychaudhuri equation, to obtain the trace of the second-order metric tensor. (Actually, in 
order to obtain the sub-leading term we also need to use the energy constraint). The resulting 
expression for (f)^ is 

^ = {-f^M + (VV) 2 ) + ^ + ^VV) • (51) 

(2) (2) 

The expression for Xslj ls obtained by first replacing (jy s ' into the remaining equations and 
solving first the energy constraint, next the momentum constraint and finally the (traceless 
part of the) evolution equation. We obtain 

(2) T 4 „ ^ _„ o . , NO _ 19 



Xsij 12g 
5r 2 



19 \ 



+ ~y y-6<P,i<Pj - + 2ip' k ip :k 5ij + -<pV 2 <p6ij) + TTij . (52) 

The transverse and traceless contribution ir^, which represents the second-order tensor mode 
Xs}/ generated by (growing-mode) scalar initial perturbations, is determined by the inhomo- 
geneous wave-equation 

<• + \< 3 ~ V% = -^V 2 ^- , (53) 

with 

Sij = V 2 ^ % + *o,*j + 2 (v, tJ VV - ¥,ik¥' k j) , (54) 

where 

V 2 * = ~\ ((VV) 2 - w*) ■ (55) 
This equation can be solved by the Green's method; we obtain for that 

r 4 4r 2 

*ij ( x - r ) = ^ <% ( x ) + — % ( x ) + *ii ( x ' r ) > (56) 

where V 2 7^- = <Sjj and the remaining piece tt^, containing a term that is constant in time and 
another one that oscillates with decreasing amplitude just like the linear tensor modes, can be 
written as 

7f y (x, r) = j±pj rf 3 kexp(,k • x)i^(k) (I - , (57) 

with Sij(k) = J <i 3 xexp(— ik • x)«Sy(x). 



4 Perturbative approach in the Poisson gauge 



The Poisson gauge, recently discussed by Bertschinger |13[, is uniquely defined by ujpi^' 1 = 
XPij^'^ = 0. This gauge generalizes the so-called longitudinal gauge to include vector and 

(r) (r) 

tensor modes. The latter gauge, in which uj Li = x L lj = 0, has been widely used in the literature 



to investigate the linear evolution of scalar perturbations |L4] . Since vector and tensor modes are 
set to zero by hand, the longitudinal gauge cannot be used to study perturbations beyond the 
linear regime, because in the nonlinear case vector and tensor modes are dynamically generated 
by the growth of scalar modes. 

Instead of writing the perturbed Einstein's equations directly in this gauge, we will trans- 
form the synchronous gauge quantities, by means of the second-order gauge transformations 
introduced in This problem has been dealt with in detail in Refs. [T], [12]]. 

4.1 Transforming from the synchronous to the Poisson gauge 

In these notes, we are only interested in calculating first and second-order perturbations of the 
metric tensor. A gauge transformation of order r can be associated to a coordinate transfor- 
mation along r independent vector fields 0^, given by 

= + \ - %) + " " " ■ (58) 

The general gauge transformation rule for a tensor T, with background value To and fluctuation 
AT = £~ i(l/H)5T is 

5f = 5T+£ kl) T , (59) 

to first order, 

5 2 f = 5 2 T + 2£ kl) 5T + T + £ i{2) T , (60) 

to second order, etc.. Here £^ is the first-order Lie derivative along the vector £| is the 
second-order one, etc., and we remind that the Lie derivative acting on a covariant tensor T^ u 
of rank two (which is the case we are interested in here) is 

(61) 

A general r-th order gauge transformation is then determined by a set of r vectors £Y r \, whose 
components we can also split in scalar and vector modes, as usual 

£?,) = « W , (62) 

and 

£* r) = d'p^ + d {r)i , (63) 

with didW = 0. 

We are here only interested in performing the gauge tranformation on the covariant metric 
tensor a 2 {T)g^ v . In transforming from the synchronous to the Poisson gauge, to first order, one 
obtains 

V£> =a«' + -aW, (64) 

T 

a« = , (65) 
<*?\ = , (66) 
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D y ( X f )» + 2/3«) = , 

Y (l)-L , rf « _0 
X s + a {i,j) - u > 

Ap AS ' 



(1) 



(67) 
(68) 
(69) 
(70) 



From these equations, given the first order perturbed metric in the synchronous gauge, we can 
obtain the parameters of the gauge transformation and the perturbed metric in the Poisson 
gauge. 

The second-order transformation, in the particular case where only scalar modes are present 
a)±_(i)T_; _ n . — , . ,(i) (i)± 



at first order (x s [j = X s ij = v\ 1)± = 0, implying d\ = uy t = x P \j = 0), yields 



(2) 



4 2) = # + /3(d [2 + ^?>) - - 

- 1 (-40« + /3; 1} 9 + p( x) di + + ^V 2 /3 (1) ) V 2 /3 (1) 

«.(!) 2 «(i)A , 2 ,3(i)^i 2 



+ /% l 2 C - + - -«(2) - gV^) , 



(71) 
(72) 



(73) 



X% = X% + 2 (Jv'Ad - 40« - f3' {1) d - D y /J (1) 

- 4 toft, - k/^g) + 2 (45) + D ^ (2) ) • ( 74 ) 



One can, at least implicitly, compute the parameters involved in this second-order transfor- 
mation in terms of the second order perturbed metric in the synchronous gauge and first order 
perturbed quantities, namely 

V 2 V 2 /3 (2) = -\ X fff + 6tf>*7$ - 2V 2 0«V 2 Ai) + 80pV 2 /^ } + 40«V 2 V 2 /3 (1) 



+ 4V 2 /$/$ - - V 2 /%V 2 /#> + ^/^ - =V^ (1) V^V^ (1) + 



2/3 V72V72, 



^V 2 /^ V 2 /3 ( ' 1} + 2/^V^' + faW^i) + %^Vff' , 



frfoW' , ^2^2 ,o/ 



72^2^(1) 



(75) 
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(76) 



V 2 a (2) = V 2 /3; 2) -2(20p + ^; + iv 2 /3j ) )/3f 



2 20« + fifa - -V 2 /3 (1) V 2 /?^ - 2/3g/3; 



0(1)' 



(77) 



4.2 Poisson gauge results 

Using the metric perturbations in the synchronous gauge presented in Section 2, we can now 
fix the parameters of the first-order gauge transformation as 

a« = ^, (78) 

= (79) 

and S 1 ^ = 0, because of the absence of linear vector modes. We can then compute the first- 
order metric perturbations in the Poisson gauge, 

ij? = <$ = ip, (80) 
Xpij - Xij ■ (81) 

These equations show the well-known result for scalar perturbations in the longitudinal gauge 
and the gauge invariance for tensor modes at the linear level. 



The parameters of the second-order transformation read 15 



where 



and 



« (2) = -^ + r(H^ + 4e ) , (82) 
ff» = ^(^V,-^o)4(^ + 6B ), (83) 



V 2 6 = * - W<P,i (84) 



V 2 rff = r 2 (-^jVV + - f*«w) • (85) 



For the second-order metric perturbations, one obtains [IS 

# = r 2 (i^V, - ^o) + fv 2 + 129 , (86) 
A 2) = r 2 Q^'V, - ^o) + \p 2 ~ 86 , (87) 

V 2 4 2)i = -|t (^VV - ^>,; + 2*j) , (88) 
X% = Try. (89) 
Quite interesting is to write down the equation governing the evolution of the second-order 

(2) 

tensor modes, Xplj — ^iji m the Poisson gauge. It reads 

t'i + ^ij ~ V% = ~%i , (90) 
with the source %j defined in Section 3.2 [after Eq. (56)]. 

(2) (2) 

Note also that the expressions for tjj p and (j) P could have been recovered, except for the 
sub-leading time-independent terms, by taking the weak-field limit of Einstein's theory (e.g. 
ref. |[L^| ) and then expanding in powers of the perturbation amplitude. 



5 Discussion 



As we have seen, both in the synchronous and in the Poisson gauge, second-order tensor modes 
are dynamically generated by the gravitational instability of scalar fluctuations. The form of 
these modes is however quite different in the two gauges. The synchronous-gauge tensor modes 
contain four terms: the first one, oc r 4 , can be easily seen to represent a Newtonian contribution, 
describing the dynamical tidal induction acting from the environment on the fluid element, then 
there is a post-Newtonian term, oc r 2 , a constant post-post-Newtonian term, required by the 
vanishing initial conditions, having no obvious observational effects, and, finally, a wave-like 
piece, which has just the usual form as free cosmological gravitational waves. Quite interesting 
is the fact that the Newtonian and post-Newtonian terms are dropped by the transformation 
leading to the Poisson gauge, so making the physical meaning of our second-order tensor modes 
more transparent. 

It has been argued that a possible observational evidence for these tensor modes could 



be in terms of a tensor generalization of the so-called Rees-Sciama effect [17, 18, 19, 20, 21 



23, EM|, a secondary anisotropy of the Cosmic Microwave Background. Let us discuss this 



point in more detail. Mollerach & Matarrese [15| have recently obtained a general formula for 



the full (scalar, vector and tensor) Rees-Sciama effect. Their expression reads 

1 r Xe / , mv ,mv OV a 1 



ST RS = - I d\ 



2Jx 



(ViW' + ^+^V-ix^eV) , (91) 



where e* represents the unit vector of the unperturbed incoming photon and the integration is 
along the unperturbed photon path. The first two terms correspond to the scalar contribution 
to the effect, while the third and fourth ones to the vector and tensor contributions, respectively. 
This splitting of the contributions depends strongly on the gauge choice. For example, it can 
be shown that in the Poisson gauge the contributions coming from the vector and tensor pieces 
are subdominant with respect to the scalar ones, while in the synchronous gauge the scalar and 
tensor contributions are of the same order of magnitude. 

One may then wonder whether there is any hope to detect the cosmological stochastic 
gravitational-wave background produced by second-order scalar fluctuations, which we have 
discussed in these notes. The general prospects for detecting a stochastic gravitational-wave 



background have been nicely reviewed by Bruce Allen at this meeting (see also Ref. It 
is, of course, the oscillating part of xij which is of relevance for earth or space detectors. The 
problem for these wave-like modes is that their energy density suffers the usual a -4 dilution 
caused by free-streaming inside the Hubble radius, while on the horizon scale their closure 
density is already extremely small, Q gw ~ 5% (where Sh is the rms density contrast at horizon- 
crossing), because of their secondary origin. More promising is the possibility that a non- 
negligible amount of gravitational radiation is produced during the strongly non-linear stages 
of the collapse of cosmic proto-structures. This issue would however require a non-perturbative 
approach, which is beyond the available analytical techniques in this field. A first step in the 
direction of estimating the amplitude of such an effect has been recently made by Matarrese 
& Terranova ||, who studied the collapse of homogeneous triaxial ellipsoids embedded in an 
Einstein-de Sitter universe, including post-Newtonian tensor modes, and showed that these 
post-Newtonian terms tend to dominate over the Newtonian ones, during the late stages of 
collapse. 
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